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ABSTRACT 

We perform a blind search for the variability of the 7 -ray sky in the energy range 
E > 1 GeV using 308 weeks of the Fermi-LAT data. We use the technique based on 
the comparison of the weekly photon counts and exposures in sky pixels by means 
of the Kolmogorov-Smirnov test. We consider the flux variations in the region signifi¬ 
cant if statistical probability of uniformity is less than 4 x 10“®, which corresponds to 
0.05 false detections in the whole set of 12288 pixels. Close inspection of the detected 
variable regions result in identification of 8 sources without previous known variabil¬ 
ity. Two of them are included in the second Fermi LAT source catalogue (FBQS 
J122424.1-f243623 and GB 6 J0043-f3426) and one (3EG J1424-f3734) was reported 
by EGRET and also was included in the First Fermi LAT source catalogue (IFGL), 
but is missing in the 2FGL. Possible identifications of five other sources are obtained 
using NED and SIMBAD databases (IRXS J161939.9-b765515, PMN J2320-6447, PKS 
0226-559, PKS J0030-0211, PMN J0225-2603). These new variable gamma-ray sources 
demonstrate recurring flaring activity with time scale ^ weeks and have hard spec¬ 
tra. Their spectral energy distributions deviate significantly from a simple power-law 
shape and often peak around ~GeV. These properties of activity are typical for flaring 
blazars. 

Key words: methods: data analysis, gamma-rays: general, BL Lacertae objects: 
general, galaxies: active 


1 INTRODUCTION 

In the previous paper llPshirkov fc Rubtsovll20l^ ) (Paper I) 
we developed a method for the blind search of long time- 
scale variability in the Fermi-LAT data. We binned the 
gamma-ray sky into a large number of relatively small pixels 
(~ 3.6 sq.deg.) and, using Kolmogorov-Smirnov test, com¬ 
pared a cumulative distribution functions of the registered 
photons and exposure. If there is no significant variability in 
a pixel these distributions should be close. The statistically 
significant deviation of the flux from uniform distribution is 
considered as an indication of the presence of some variable 
source in the corresponding pixel. The method is sensitive to 
the slight gradual changes in the flux that are more difficult 
to discover with the alternative techniques. However, that 
comes at a price of somehow reduced sensitivity to short, 
burst-like events. 

In the Paper I 168 weeks of data were analyzed, 151 


* E-mail: pshirkov@sai.msu.ru 
f E-mail: grisha@ms2.inr.ac.ru 


pixels with a significant variability were found (see Section 
[2]for the definition of significance). Variability in these pixels 
was due to 118 sources (some sources contributed to several 
neighboring pixels), 27 out of which had not been previously 
reported as variable in the Astronomer’s Telegram^( ATels), 
on the Fermisky blofj3 or in the dedicated papers. 

In this paper we extend the data set by 140 weeks to the 
first 308 weeks of the Fermi mission livetime. This reanal¬ 
ysis is justified not only because of almost twofold increase 
in the statistics but also by using the new reconstruction 
framework Pass 7 Reprocessed ( P7REP) recently pres ented 
by the Fermi-LAT collaboration jBregeon et ahlBoiSl ). The 
new improved instrument response functions (IRFs) may 
result in a change of calculated exposures and thus signifi¬ 
cantly affect the results of our analysis. 


^ http://www.astronomerstelegram.org/ 
^ http://fermisky.blogspot.com/ 
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2 DATA AND METHOD 

In this work we use the LAT Pass 7 Reprocessed weekly all¬ 
sky data publicly available at the Fermi mission websitj^- 
The analysis spans the time period of 308 weeks from Au¬ 
gust 04, 2008 to June 23, 2014, corresponding to the mission 
elapsed time (MET) from 239557417 s to 425179643 s. We 
use the ’SOURCE’ event class photons with E > 1 GeV 
and impose an Earth relative zenith angle cut of 100° and 
rocking angle cut of 52°. The data selection criteria and our 
method are identical to the ones described in details in Pa¬ 
per I, here we only briefly repeat the main points: 

• We bin weekl y data files using HEALPIX package 
(|G6rski et al.l[20n5h into a map of resolution Wide = 32 in 
galactic coordinates with ’RING’ pixel ordering. Total num¬ 
ber of pixels is equal to 12288 and the area of each pixel is 
3.6 sq. deg, chosen according to the size of Fermi-LAT point- 
spread function (PSF) above 1 GeV which is approximately 
1 °. 

• We estimate an integral weekly exposure for each pixel 
using the standard Fermi-LAT tools gtltcube and gtexpcube 
from the Science Tools v9r32p5 package. 

For each pixel we count the number of photons in each of 
308 weeks and consider corresponding values of weekly ex¬ 
posure. Cumulative distribution functions (CDFs) 'P{t),£{t) 
for both photon counts and exposure are constructed. 

• We compare distributions in each pixel using the 
Kolmogorov-Smirnov test and find all pixels with a prob¬ 
ability Pks lower than the threshold value Po = 4 x 10“®. 
This threshold leads to an average 0.05 false detections in 
the blind search through the whole 12288 pixel set. 

We have produced a dedicated set of Monte-Carlo sim¬ 
ulated photons with the gtobssim utility using the positions 
and s pectra of sources from the 2FGL catalogue l|Nolan et al.l 
l2012h and the latest templates for the diffuse Galactic and 
isotropic gamma-ray radiation. All sources were naturally 
simulated as non-variable. Repeating the analysis on events 
from this Monte-Carlo set we arrive at the lowest Pks value 
around 10“"* in total agreement with the expectations. 


3 RESULTS 

The analysis resulted in 231 pixels with Pks < Pq. We com¬ 
pare our findings with the results of the Paper I and con¬ 
centrate mainly on the differences. First, there are 20 pixels 
reported in the Paper I, but not detected now. These changes 
could be due to two main causes: 1) variability may decay 
on longer timescales and thus the corresponding probabil¬ 
ity for this timescale would increase and, eventually, would 
overshoot the chosen threshold 2) rehned exposure calcula¬ 
tion with the improved IRFs may affect an apparent pixel 
variability obtained with the older IRFs. The most intrigu¬ 
ing case was the one of the pixel 5637 which contains the 
Geminga pulsar. We have reported Pks = 2.3 x 10“^ for it 
in the Paper I. However in our new analysis this probabil¬ 
ity was only 9 x 10“®. This decrease is mostly attributed 
to the change of the reconstruction - the probability for 


® http://fermi.gsfc.nasa.gov/ssc/data/access/ 


the first 178 weeks recalculated using P7REP gives only 
Pks = 4.5 x 10"^ 

Second, we deal with the pixels that have not been re¬ 
ported in the Paper I (100 pixels total). In order to crudely 
locate a source within the pixel to check whether its vari¬ 
ability was already known, we plot the photons belonging 
to the weeks with the highest flux and find the center of 
their clustering. After that, the attribution was hrst made 
using 2FGL catalogue and then using SIMBAE0 or neeQ 
databases if no counterpart was found in the 2FGL cata¬ 
logue. 

Finally we focus on the pixels where no sources with pre¬ 
viously known variability was found. That comprises both 
pixels with known gamma-ray sources that were supposed 
to be non-variable and pixels with some previously unknown 
gamma-ray sources not included in the 2FGL catalogue. 
Therefore we discard the pixels containing sources with vari¬ 
ability repo rted in ATels, on the F ermi sky blog, or in th e 
papers e.g. dAckermann et al]l2013h and dAbdo et al.ll201ol L 
or, finally, if their variability indices in the 2FGL catalogue 
exceed the ’ofRcial’ threshold value 41.6. In the end it left us 
with 8 pixels out of 100 corresponding to 8 different sources. 

With 8 remaining pixels we perform more thorough 
search in the following manner: first, we find a location of 
the apparent excess using the gttsmap utility. Known 2FGL 
sources and the Galactic and isotropic diffuse components 
were included in the initial sky model that was provided as 
an input for gttsmap. In two cases (pix.34, pix. 9067) no new 
significant excess was found meaning that the flares occurred 
in the sources from the 2FGL catalogue. These sources were 
identified by their large TS values in the results of stan¬ 
dard gtlike ht. Variability in these pixels was caused by 
2FGL sources FBQS J122424.1-b243623 (MS 1221.8-f2452) 
and GB6 J0043-I-3426 correspondingly. The hrst one is BL 
Lac type object at the redshift z = 0.218 that was recently 
observed al so in the TeV energy range by the MAGIC col¬ 
laboration (ICortinall2013ti . the latter one is a blazar at the 
redshift z = 0.966. In the remaining six cases (see Fig. [TJ 
we have used gtfindsrc in order to rehne the positions of the 
sources (see Table [T]). 

We have performed searches in the NED and SIM- 
BAD databases using the best ht positions and reported 
errors obtained earlier. One of the sources may be identihed 
with a high conhdence: the reported variability in the 431- 
th pixel was caused by the 3EG J1424-I-3734 BL Lac type 
object that wa s included into the 3rd EGRET catalogue 
llHartman| []lI9£ ) and also into the hrst Fermi LAT catalogue 
(IFGL) I Abdel 2010l j. However it is absent in the 2FGL cat¬ 
alogue and that clearly indicates its long-time variability. 
The other hve sources were linked with the known radio- 
and X-ray sources from the NED and SIMBAD databases. 
One of the hve, PK S 0226-559 (pi x.11255), is a blazar with a 
redshift 2 = 2.464 (iHealev et al.l[200^1 . Two other sources, 
PKS J0030-0211 (pix. 11696) and PMN J0225-2603 (pix. 
119010) are hat-spectrum radio sources. Two remaining can- 


^ http://simbad.u-strasbg.fr/simbad/ 

® http://ned.ipac.caltech.edu/ 

® There is a possibility that the variability in this pixel was caused 
by another source, IRXS J022517.5-260926 located only 6.6 ar- 
cmin away from our best localization of the source but a flat- 
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didate sources IRXS J161939.9+765515 (pix. 2663) and 
PMN J2320-6447 (pix. 10871) has no classification type. 

We have accomplished a detailed analysis of the spectra 
of these sources. We have carried out gtlike fit in two energy 
intervals - (0.1-100 GeV) and (1.0-100 GeV) adding new 
sources to the source model files obtained from the 2FGL 
catalogue. Also we conduct fits in (0.1-100 GeV) interval 
using the spectral models other than just a simple power- 
law. The results are presented in the Tables I2I3I We plot 
the spectral energy distributions (SEDs) using the python- 
based package likeSED.pj[3|(figs. (12131) 1. The spectral shapes 
demonstrate significant deviation from the simple power law. 
For a majority of sources their SEDs peak at around GeV. In 
all cases the spectra are rather hard -- that makes them look 
very sim ilar to the SEDs of t he brightest flares of powerful 
blazars JPacciani et al ][2^. Thus the observed variabil¬ 
ity seems to be caused by some sources with temporal and 
spectral properties close to highly-variable blazars. 


4 CONCLUSIONS 

We have performed a full sky blind search for regions with 
variable flux at energies above 1 GeV in the 308 weeks of 
the Fermi-LAT data. The variability of 8 gamma-ray sources 
was identified for the first time. Two of the sources are in¬ 
cluded in the second Fermi LAT source catalogue (2FGL), 
one was reported by the EGRET and in the First Fermi 
LAT source catalogue (IFGL), while the gamma-ray radi¬ 
ation from the five other sources is detected for the first 
time. The latter, using the NED and SIMBAD databases, 
were linked to the sources known in radio and X-ray bands. 
All detected variable sources demonstrate recurring flaring 
activity at time scales ~ weeks which is typical for blazars. 
The SEDs of the sources in their high states are hard and 
deviate significantly from the simple power-law. The spectra 
of the most of the sources peak around few GeV. We note 
finally, that the variability property was for the first time 
employed to detect new gamma-ray sources. 
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spectrum radio source PMN J0225-2603 looks a more plausible 
candidate. 

http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 


no 

Pixel no. 

r 

b° 

error, ' 

1 

431 

67.33 

68.10 

3.0 

2 

2663 

110.32 

34.47 

6.6 

3 

10871 

317.64 

-49.67 

5.5 

4 

11255 

278.42 

-56.52 

4.0 

5 

11696 

110.79 

-64.59 

3.0 

6 

11901 

215.54 

-69.14 

3.5 


Table 1. Positions of the gamma-ray sources not included in the 
2FGL that were obtained with the gtfindsrc utility. 

UPDATE 

When this paper was already finished, the Fermi LAT T hird 
Source Catalog (iThe Fermi-LAT Collaboration! l2015l l has 
appeared on the arXiv. We have checked for our variable 
sources from the Table [4] there. Both 2FGL sources (pix.34, 
pix.9067) are in the 3FGL catalogue as well. Sources that 
were responsible for the variability in the pixels 11255, 
11696, 11901 (PKS 0226-559,PKS B0027-024, PMN J0225- 
2603) are firmly estabisished in the new catalogue. How¬ 
ever, 3EG J1424+3734 (pix.431) is absent from the cat¬ 
alogue. Also our identification in the pix. 2663 (IRXS 
J161939.9+765515) is 1.3° away from the closest 3FGL 
source, which in turn lacks identification. Our identifica¬ 
tion of the source in the pix. 10871 PMN J2320-6447 dif¬ 
fers from the 3FGL - PMN J2321-6438 (I = 317.6567°, 6 = 
—49.8059°) which is 26.1' away from our best localization 
position {I = 317.64°, 6 = —49.67°) that has an error 5.5'. 
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no 

Pixel no. 

Floo, 10 ® 
ph cm“^ s“^ 

^FlOO; 10 ® 
ph cm~^ s“^ 

Tioo 

i5rioo 

Fiooo, 10 ^ 
ph cm“^ s“^ 

(^Fiooo, 10 ^ 
ph cm“^ s“^ 

Tidoo 

iSriooo 

1 

34 

3.7 

0.04 

1.71 

0.004 

0.79 

0.04 

2.00 

0.06 

2 

431 

3.5 

0.03 

1.710 

0.005 

0.85 

0.05 

2.07 

0.04 

3 

2663 

0.68 

0.03 

1.72 

0.01 

0.14 

0.05 

2.18 

0.34 

4 

9067 

5.3 

0.07 

1.68 

0.006 

1.2 

0.06 

1.97 

0.03 

5 

10871 

2.3 

0.07 

1.64 

0.01 

0.75 

0.19 

2.44 

0.30 

6 

11255 

8.5 

0.2 

1.99 

0.007 

0.95 

0.19 

2.40 

0.29 

7 

11696 

22 

2.9 

2.11 

0.098 

0.75 

0.19 

2.44 

0.30 

8 

11901 

12 

0.1 

2.15 

0.003 

0.95 

0.16 

2.81 

0.26 


Table 2. Spectral properties of the sources in their high states. -F 10 O 5 ^lOOO ^-re photon fluxes at energies higher than 100(1000) MeV, 
(^-^100? ^-^1000 are respective errors on these fluxes, Fioo, Fiooo are spectral indices in the simple power-law model at these energies, and 
5Fioo? ^Fiooo are their errors. 


no 

Pixel no. 

PL 

BPL 

LP 

PLEC 

1 

34 

134.8 

137.6 

136.5 

136.5 

2 

431 

159.0 

162.9 

165.7 

167.5 

3 

2663 

25.0 

25.4 

26.2 

26.5 

4 

9067 

260.1 

264.9 

262.7 

261.0 

5 

10871 

43.2 

52.6 

58.6 

56.9 

6 

11255 

244.6 

248.4 

249.1 

250.2 

7 

11696 

297.9 

298.8 

298.5 

299.1 

8 

11901 

203.2 

208.4 

208.3 

210.0 


Table 3. Goodness-of-fit {TS obtained during gtlike fitting) with different spectral models: simple power-law (PL), broken power-law 
(BPL), log-parabola (LP), and power-law with exponential cut-off(PLEC). The best model for each source is highlighted in bold. 


no 

Pixel no. 

1° 

b° 

F^phot 

Pks 

source identification 

distance, ' 

source type 

In 2FGL 

1 

34 

236.25 

84.15 

610 

1.2 X lO”^'^ 

FBQS J122424.1-I-243623 

- 

TeV BLL 

yes 

2 

431 

67.33 

68.10 

419 

2.2 X 10“® 

3EG J1424-I-3734 

2 

BLL 

no 

3 

2663 

110.32 

34.47 

574 

8.5 X 10 “ 

IRXS J161939.9-I-765515 

10 

? 

no 

4 

9067 

120.94 

-28.63 

840 

7.9 X 10-12 

GB6 J0043-I-3426 

- 

blazar 

yes 

5 

10871 

317.64 

-49.670 

351 

1.6 X 10-1’ 

PMN J2320-6447 

0.2 

? 

no 

6 

11255 

278.42 

-56.52 

55 

5.9 X 10-21 

PKS 0226-559 

7 

blazar 

no 

7 

11696 

110.79 

-64.59 

326 

1.3 X 10-® 

PKS J0030-0211 

2 

fiat-spectrum radio source 

no 

8 

11901 

215.54 

-69.14 

289 

3.5 X 10-® 

PMN J0225-2603 

7 

fiat-spectrum radio source 

no 


Table 4. List of the sources with previously unreported variability. The coordinates of the sources are taken either from the 2FGL or 
from the results of the gtfindsrc searches. Np^ot is the number of photons observed in the pixel during 308 weeks, Pks fLe probability 
that the pixel is non-variable, source identifications were made using the NED and SIMBAD databases. 
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Figure 1. TSMaps for the pixels listed in the Table|4] (a) Pixel 431: TSmasj=115 (weeks 66,293,301,302,316); (b) Pixel 2663: TSmax=17 
(weeks 109, 245, 250, 260, 261, 306); (c) Pixel 10871: TSmax=64: (weeks 142, 170, 311, 314, 315); (d) Pixel 11255: TSmax=126 (weeks 
139, 238,276, 283, 285); (e) Pixel 11696: TSmax=112 (weeks 200,201,276); (/) Pixel 11901: TSmax=102 (weeks 134,135,138,184,289) 
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Figure 2. Spectra for sources from pixels (a) 34, (i») 431, (c) 2663, (d) 9067. Left: TS value for each bin of the spectrum. Right: Spectrum 
of the source in its high state. 



















































































7 


Revisiting variable gamma-ray sky at 1 GeV ... 


(a) 



(b) 


(c) 


(d) 



Energy (GeVj 



Energy (QeV) 



Energy (GeV) 



Energy (GeV) 


Figure 3. Spectra for sources from pixels (a) 10871, (b) 11255, (c) 11696, {d) 11901. Left TS value for each bin of the spectrum. Right 
Spectrum of the source in its high state. 




















































































